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■ Iabstract .... J j; ; 1 1 . . 

r S«fExacl image method, recently introduced for the exact solution of electromagnetic field 
problems involving homogeneous half spaces and microstrip-like geometries, is developed 
here for the problenl of homogeneous slab of dielectric and/or magnetic material in free i 
space. Expressions for image sources, creating the exact reflected and transmitted fields; 
are given and their numerical evaluation is demonstrated. Nonradiating modes, guided by" 
the slab and responsible for the loss of convergence of the omage functions, are considered 
and extracted. The theory allows, for example, an analysis of finite ground planes in 
microstrip antenna structures. 
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INTRODUCTION 

Exact imago theory, originally introduced for the Sonimerfeld half-space problem fjl, 
;2' was later generalized to the microstrip geometry [3], |4j consisting of a dielectric slab 
plus a perfectly conducting ground plane, To account for the more realizable case of finite 
ground plane, it is necessary to further extend the theory to the dielectric slab geometry, 
because the currents in the ground plane can then be solved through integral equation 
formulation. This problem was recently worked upon by Mosig |5] by using a numerical 
method to evaluate the ill-behaved Sonimerfeld integrals, Because such integrals could be 
avoided in earlier applications of the exact image theory, the same principle suggests itself 
also for the present geometry, 

The present paper is a continuation to the previous microstrip paper [4] with a minimal 
change in notation. The medium in the present case, however, is assumed to be more 
general, with both electric and magnetic susceptibility, and also the image giving the exact 
field transmitted through the slab will be developed. Because of this, the expressions have 
grown in extension and complexity in comparison to those of [4], The complexity in image 
expressions is not, however, a major problem in the field calculation, because the image 
functions need only be calculated once for storage in the computer memory, after which 
their repetitive use in solving integral equations in the presence of the slab is possible with 
the present method. 


BASIC THEORY 

The solution of the dielectric slab problem can be obtained by decomposing thg original 
source into symmetric and antisymmetric parts, which correspond to problems involving 
microstrip- type geometry. In fact, because of the symmetries of the corresponding fields, a 
plane of perfect conductor or magnetic conductor can be placed at the plane of symmetry 
in the middle of the slab and thus the problem is split up in two parts, one with the 
microstrip geometry and the other with the magnetic conductor microstrip geometry. 

It is, however, preferable to make a fresh start for the slab geometry, because it is 
no more' complicated than the study of magnetic conductor microstrip problem. Also, we 
arc going to study the general dipole source above a slab with both electric and magnetic 
discontinuity and not just a transverse dipole at a dielectric interface as was done in [4]. 
However, the start with Fourier transforms similar to that in [I, part III], [4] will not 
be repeated here. The notation w’ill be compatible with that in }4j, with minor obvious 
changes. (( 

The reflection dyadic for the electric field transverse to the z direction can be easily 
written from the partial wave reflection geometry, Fig. 1, from the knowledge that. TE 
and TM polarizations are its eigenvectors; 


S = := (1 - • 
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with f( — I — nil. the two-dimensional unit dyadic and r is the interface reflection dyadic. 
The inverse dyadic must be understood as a two-dimensional inverse [8]. We can further 
write 
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Here, A' is the two-d imensio nal Fourier variable on the plane transverse to the u vector 

and /? =i s/k- ~ A' 2 ,/3j = \Zfiek 2 - K 2 .are the propagation coefficients of the Fourier 
transformed wave in z coordinate; (direction ii).: Note that the thickness of the slab is 
denoted by 2d to be compatible with the tnicrqstj’ip substrate thickness d of [4], 

Also, the field transmitted through the slqb can be calculated from the knowledge of 
the total transmission dyadic, which is of the form 

f = r ™*||»« +r ™« (5) 

t TE>TM _ £ ^ r TE,TM )2n _ ( r TE,TM fn+'l\ e ~jfa(2n+l)2d ^ 

n=0 

From (2), (6) there results a relation between the reflection and transmission coeffi- 
cients, and (5) can be given different dyadic forms: 

f = (7/ - 1 • ?) e -tf' 2d = e~iP' 2d (7 i - F 2 ) . (7 t - re~tf' Ad y l = i(I,„ - %), (7) 

where R m , R e are the reflection dyadics of the microstrip geometry with magnetic and 
electric conductor ground planes^ respectively. As a check we see that for F 2 — 0, or = It, 

W& have respectively from (7), T — Itc~^ l2d , or — 0, which correspond to the cases of 
perfectly transparent and perfectly reflecting slab. 

The* reflected and transmitted transversejelectric fields can be written in the form given 
in [1, p.1028] for the general electric dipole J(r) vIL6(r — uh) as follow’s: 

_ h — € -3fH^h) 

e(A’,5) = -(a + jpb)-R — 2 j — -, for ; > 0 (8) 

- =e jf3{i~h+2d) 

c{K,z) ~ -(a + j/J6) • T ■ , for z < ~2d (9) 


a = jk 1 r)(k 2 It - KK) - vIL, b — k~*T)KJL{u • v). (10) 

To calculate the image for a transversal source, the reflection dyadic can be split up 
into two parts as was made in [4]: _ 

1 = R TE 7t + Ro 1 ^-, (11) 
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where we have 


Ro ~ %(R ni 
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To be able to construct image functions, the reflection cocfficents are first expressed in 
terms of the reflection variable r, (which equals r TE for // = ! and for e = 1): 


/?+/?r 

(15) 

Thus, we have 


Ti? A/ + r 
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(16) 

TM E + r 

1 + Et' 

(17) 

il 

vt 4/1 1 \ 

0 p — c \ 1 4* A/r 1 4- Er) ‘ 

Cis) 

The reflection coefficients W ere given as power series of , 

For Ro the following series can be written from (12): 

r TAf in (2). 

OO 00 

^0 = (1- c-^ 4rf )r 0 £ jr Kmn(r TE ) m (r™ fc^m+n^ 

(19) 


m —0 ti=0 


where K m n equals 1 when both m and n are even or odd, and = 0 otherwise, or explicitly, 
K m „ = jl — |-1 ) ,n ‘ tn ;/2. Another form, less symmetric than (19), but better suited to 
successive computation of terms, is the following: 

Rq - ro ‘ - c "i/ ? i(?+l)4<fj _ 

q=0 i=0 r— 

Whichever form of (19), ( 19a) is used, (16) and (17) can be substituted to obtain repre- 
sentations in terms of powers of the reflection variable r, 

Like in the previous papers [l] - |4], the reflection coefficients will be presented as 
integral transforms of Laplace type Ito obtain image source functions. To this end, the 
following integral identity, invoked from [4], will be applied: 

r m e-^ c * [ F m (c,t ) + 6 m0 6(i - c)J (20) 

JO 
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'Here and henceforth, the delta fund ion />(»), also denoted />+(,r) in il;, is defined to be 
zero for :r < 0 so that 6(x)dx — 1. F m (cJ) is the following function; 


F«(c,t)==~-I l {T)U(t~c), 


(21) 




F, n (c,i) * ~{(i- c)h m+ ,(r) - (t + c)/ 2m _,ir)) W - c) 


IB 

2 


W - c)) 2m+1 ^giTr - W - c )) 2w " , ^£rr] W - *)’ ™ > 0. 


(22) 


with 


T, r « B\/l 2 -c 2 , 


(23) 

U(t) is the Heaviside unit step function making each function satisfy F m (c , f) = 0 for /. < c, 
and / n (r) denotes the modified Bessel function. For c = 0 (21), (22) reduce to 


0 77| 

mo p-o, F m (o,o = -j-hmmm. 


Further, we need the following expansions [4] 

r _ 


Jt=0 
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*=o 

where the coefficient function is defined by 

71 |A' 
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(24) 

(25) 

(26) 

(27) 


q°(*) = %, c 0 n (*) = (-*)”> 


C’/(*) = (~x)* +I - (-*)*- 1 , C'”(x) = n[( r x) n+i - (-x) 

Here, n\k denotes the smaller of the numbers n and k. Values for the function can also 
be computed by reflection: 


k— 1 fin i 


\n-li 


C^) = -Ct(x), k,n> 0, 


(28) 


and using the recursive formula 


cm = -icf-'f.) - ^4 E<-*) i - p cr , (*>- 


p=0 


X 


(29) 


(29.) is more rapid (hail (27) if many values of Cj!(r) are to be calculated in computer 
memory for successive values of k and n with the same value of the argument x (constant 
f or //). The coefficient function Cj!(r) can also be defined through the expansion 



00 

fanh n (0 + ('}! ( Ian h0)1anh k ‘ <jt 

lixQ •' 

\\ ■ 


VED REFLECTION IMAGE 

Lei us first construct the image due to the vertical dipole J(f) — uIL6(p)6(z — h), 
which lies at the distance z — h from the slab. The transversal electric field in Fourier 
space is obtained from (8) when v ~ u: 

HK,z) = -jpk-'nKILR™ ~ , (30> 

If the reflection coefficient, is written in the integral form Jj 

• • . ' / • ' i ,\ 

POO 

R 1M ^J f™(t)e~ i0l dt, (31) 

the field can be easily seen to arise from the equivalent image source [l] 

Ii(t) ~ ILf™ (t). (32) 

To find the image function thy expansions (2), (17), (20) and (25) must be applied: 

OCr OO po 

Ji™ ~ V C'l(E)r k ' +J2Y, (cl m+1 (E) - Cl m -\E)) r k e ~ jP]4md , (33) 

t’=:0 m=l k=zO 

whence the following expression for the image function is obtained: 

f™U) = C’l(E)6(1) + f; (C% m+1 (E) - Cjf m ~*(E)) 6(1 - 4m d) 

+ £ Cl(E)F k (0,t) + f; f; (cl m+1 (E) - C!l m -'(E)) F k (4md,t). (34) 

k~\ m=U= 0 ! .. r i .. i - 

|' r ■ f 

Thus, the image function consists of a series of delta functions, corresponding to dipoles 
with interspaces 4d, which is double the slab thickness, and a continuous source function. 
The function depends on both c and /x through the variable B in the functions Ff,(x,l), 
and on e through the functions G£(E). It is seen that the dipole images depend only on 
<• Taking ( — 1 , or what is equivalent, E — 0, all the delta terms in (34) vanish because 
of the following limit of the function C£(x) as x —> 0 : 
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The vanishing of delta terms is also seen if (34) is written in t he form 


f™w :z + pTi E ( hrl ) Hl ~ 4m(l) {M(i 

m - » 

f\k(-^r)* !&£ t] VW t E E “ (T>'(E)) I<\.(4md.l) 

" 1 • • V f + 1 7 1 


To obtain this, (24) has been applied. For < — 1 , all coefficients of the delta functions are 
seen to vanish, Tb< remaining image function can be written as 


f f™{t) * ~Fi(0,0 - E (F2„, +1 (4md ) <) - F 2m _,(47m/,<)). (36 ) \ 

ji. . m - l / 

This expression can be compared with that given in [4] for the TE image with f-t = 1 in 
the inicrostrip problem. 


VED TRANSMISSION IMAGE 

The corresponding t ransmission image function can be obtained from the expressions 
(6), (20), (25) for T t ™ . Thus, writing 

rc <%> 

h™(t)c-Wdt = E ((r™) U ~{r™) 2n+2 ) c-^< 2 " + ')2rf 
n=0 

00 00 

*EE (cfW - Cf +2 (F)) r*e-^t(2«-f l)2rf, (3 7) 

11=0 A’=rO 

we have 

oc oo 

* TJW (0 = E E KW - Cl n+2 (E)) [F k ((2n + l)2d,0 + 6 k0 6(t - (2n + l)2d)] 

n=0 k—0 j • ‘ 

= rTT^E £?2n ^-( 2 " + 1 ) 2rf ) + E E faPW - Cf +2 (S)) F*((2n + l)2d,<). 

^ £+ n=0 n=0Jt=0 

/ (38) 

Similarly to the case for the reflection image function, if e — 1 the delta function 
(dipole) image terms vanish except for the 'first one, which corresponds to the original 
dipole source. The resulting function for c = l ( E — 0) is of the form 

h™(t) = Hi - 2d) x ElF 2 „((2n. + l)2d,f) - F 2n+2 ((2n + l)2d,/)j. (39) 

.V^" ' »2=:0 

7 
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I he 1 rails Verse component of (lie TM electric field transmitted into the region s 2d 
cap In' written from (0) * 

- | mu e M*~h) 

r ( A , ; ) a jfik * V Kl LT tai — — ~ , (40) 

ima < >e curm t' S as (8" I* (40) cap lie interpreted as the field arising from the 

m = iLh™a), /( ( 4 i) 

where h™(t) is of the form (38). Obviously, the transmission image has no counterpart 
m the micros* rip problem, 


II ED REFLECTION IMAGE 

The image current due to a horizontal (transverse) electric dipole 3(r) = vIL6(p)6(z ~ 
//) can be divided into a horizontal and a vertical part [1], [4]. The horizontal part, is 
obtained from the function f rE {() defined through 

^=j^,TS it)e -m it , ( 42 ) 

Again, the reflection coefficient will be written as a power series of r: 


OO 00 


rIE * “• £ <i W)r k ~ £ £ (c'l m+ \M) - (43) 

*•“0 TO=1 k-o '* 

This leads to the following expression for the image function ° 

. . . . ... OO 

/' (') 5:5 ~ £ (c^\M) - Cl m ~ x (M)) 6(i ~ 4 md) 


oc 


m=t 
oo oc 


- £ C'l(M)F h ( 0,f) - £) £ (c*| m+1 ( A/) - F k (4md,1), (44) 

m=l k—0 

or, what is equivalent, 




2m 


(44a) 


8p 


£ £ (Cf'" +, (A/) - C, 2 ™- 1 !*/)) F t <W,(). 

*-l ni— 1 t-0 
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This is very similar to the function P ^ (t) given above in (34 ), In fact, changing < and 
ft, ami the sign of the function the function P ‘^ (>) is obtained, and vice versa, 

Thus, the same procedure of calculation can he used for both functions, which simplifies 
the use of the method, Tor //••- I with (85), the expression (44) reduces to 

«■ 

P' L '(n - ^((U) + Fi,„. )[4ntil,n], (45) 

wi- 1 

which is in concord with that given in ‘[4] for the microstrip geometry. In fact, a comparison 
shows us that (45) contains every second term of the inicrostrip image function, namely 
those with odd indices in the functions /^(c,t). For a microstrip with magnetic conductor 
ground plane, the sign of every second term (the even-indexed function terms) in the image 
function series is opposite to that of the normal microstrip, Thus, (45) is the average of 
the corresponding functions of the two microstrip problems. The microstrip image delta 
term due. to the ground plane is absent in the slab geometry, because it is cancelled by the 
magnetic conductor term. 

For the vertical part of the image current, a more complicated function /off) is needed: 

n m<r 0t dt> (46) 

JO 

Applying (20) in (19) with (25), (26) inserted and defiling 

+ £«(*)= 0, (47) 

the following expression for the reflection coefficient Rq is obtained: 

* = + 

Xp+ f _ c ~j/3 J (m+n+2)2<i'j _ ( 4 gj 

The symmetric double sums can be calculated more practically in the following way: 

gsO .*=0 

Thus, the partial image function /o(<) can be written in the form 

OC OC OG OG / 2 1 2 \ 

md = E E E D-u^mn (S—muietim + ~—,efmme)) 

771=0 n=0 A~0 /=0 ^ M ' 

X [Fk+ld™ + n)2d,t) - F*. + /((n? -f n + 2)2d,f)| , (49) 

or, corresponding to (48a), in the form 
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/«<')■« E E - mi - DHoj 

9=0 »=0 


■ ft (£^ D T;im'p ) + • («») 

The complicated appearance of the function /o(t) is counterbalanced by the fact that the 
functions F a (c,1) are zero for l < c, whence only a few terms must be actually calculated 
for low t values and, in many cases, the calculation can be stopped before the number of 
terms grows very large. I--.-, 

For the special case p — 1, or M — 0 the above function is simplified to 


oo oo ot> 


/o(0 = -7- E £ E KmnD?(E) [F m+ i{(m + n)2d,t) - F m+l ((m + n + 2)2d,t)] 


n»=0 n=0 /~o 


= 1 ~ L E t) - F$((q + l)4d,<)] £ D\+i^ q (E). (50) 

9=0 *=0 i =0 

In (50), only terms with s + i — 2 q > 0 contribute. This expression is in accord with the 
corresponding result (function <?(<)) for the microstrip geometry given in [4], fr 

The HED image current expression can be obtained from the combination of the hor- 
izontal and vertical partial image current functions f TE (t), /<,(/), from a comparison of 
the resulting transverse electric field, 


__ it i.' Jt-]j3(z+h) 

e(A’, z) = -jkq(It - k- 2 KK) • Ro~ • vIL - - , (51) 

with the corresponding Tree-space expression. This was done in [1, p.1029) for the Som- 
merfeld probleih and in [4] for the microstrip problem and the resulting expressions for 
the image current source are essentially the same as here, with a difference in the choice 
of int egrat ion variable in .1 [] and notations for the image functions. Thus, we can write for 
the image of nay transverse source satisfying u • J{r) == 0, 


Uf,t) = f TE (t)3c(r) - M-*fo(t)V • Jed), 


=1.-2 ft t 


IcdY- c- J(C-f), C ~ Ii - 2uu. ; (53) 

It is seen that the effective function in the vertical current component is not /o(t), but its 
derivative /^(t), which is in accord with the results of [1], [4]. 



HED TRANSMISSION IMAGE 


The expressions for the transmission image for the transverse, or horizontal, source can 
be obtained much in the same way as above by seeking the horizontal and vertical partial 
images separately. The transmission dyadic is of the form 


k m i;i™ Kh 

T « T 1E J i r T 0 ~~ , 


m 


Expressing the first transmission coefficient as 


*00 " OP 

T TE = / h^We-i^dl = £ ((r TE ) m - (r T£ ) 2 "+ 2 ) e -^i(2n+l)2rf 
^ n~G 


00 00 


£ (e$*m - cl n + 2 (M)) 

n=Ok~0 ' 

we can write for the .horizontal partial transmission image function, 


( 55 ) 


h TE (t ), 


4 ; «> n ' • 

u T ' n=0 


00 oo 


+ E E (eftM - cf +2 (»>) n«2» + 1 >2d. f )- 

n=OJt=rO ! ' 


(56) 


This is the same function as h T ^Ui) if the change e uis made, as is seen from (38). For 
fi = l the forn: ‘ ' “ 

The vertic 
coefficient Tq: 


i iuo io I'A&v same luiivuuu as u 11 uic iiiau^c t. p 

fi = 1 the form of the expression (39) is seen to result. 

The vertical part of the HED transmission image is obtained from the transmission 

7 E i rri . 


k- \ 


,—jfh 2<t _ JP. ? 


fOO 

To = / ho(l)e~ J ^ t dl ‘ 

./0 ' 

\ ( r T£')2_ (r TAi ) 2 

/ (j _ ( r rE)2 e -j/J]4rf)(i _ ( r TA;)2 f -j/3]4rf) 


Op oo 


= -ro E £ ((r™) 2 "(r T£: ) 2m+1 + (r™ ) 2 «+ J (r T£ ) 2m ) 


n=0 m=0 


- X ^ e -jA(2B+2m+l)2rf _ ( ,-'}P\ (2n+2m+3)2<f j 
Hence, the image function can be written in the following two forms 
oo oo oo oo 2 

MO = - E E E E - ^"“(JIOCf+OE)) 

«=o m=o t=d /=o ' 

2 

[c% m+l (M)Df n {E) + C'£ TO (M)D ; 2n+1 (£)) ] 


( 57 ) 


x|F* +/ ((2m + 2* + 1)2 d,l) - F* + ,((2m + 2n + 3)2</,f)j, (58) 

Qo oo 

E £[$$& + l)2rf,0 - F # ((2 9 + 3)2d, <)] 

9=0 *=o / .... o 

. • f j V 

* 1 1 [&k ££7V)<TW) •; 

«=0;=0 "' r ’ 

"^^(^T^ , <")^< , ) + '°£r*(") | f +, <«)]- . ( 58 «) 

The result has a bit more complicated appearance than the corresponding reflection 
image function fo(t) in (49), (49a). The difference arises from the difference in the partial 
wave structure^ Fig.l. For the special case p =-. 1, (58a) reduces to 


£ £[ft((2, + 1 )2d, /) - F a ({2q + 3)2d, <)) £(I>« +1 (E) ~ D?(E)}. (59) 


< 7=0 4=0 


GENERAL IMAGES 


To obtain the reflection and transmission images of the! general dipole, previous results 
can be combined. In short, for the general three-dimensional source function J(r) the 
reflection image is 


ji(r,1) = ( + f TE (t)Jt) ■ Mr) - uk- 2 f' 0 (t)V t . Mr). (60) 

The reflection field can be obtained everywhere in the half space r > 0 through an inte- 
gration of the free-space Green dyadic multiplied by the image: 


(r) = -~jki] f f G{f - f' + ut) - ji(r ' , 

i Jo Jv ~ 


If the original source is a dipole or a planar current, source of the form J(r)’ — J a (p)6(z — h), 
the integration variable t can be disposed of, because in the field integral (61) any function 
of the form G(r-r' +ut)<f>(1 )6(z' +h) can be replaced by the function —G(f )<f>(~(z' -\-h)), 
leaving out the 1 integration and taking the z' integration from -h to — oo. The image 
of the dipole consists of a set of dipoles and partially overlapping line sources along the 
negative real c axis, which start at the points +(/i -f m4d), m = 0, 1,2,... . 

For the general transmission image we have; correspondingly, 

= (h™(1)uv + h TE lA • j(f) - uk~ 2 h' 0 (t)Vt ■ J(r), (62) 

' ' (t 

which produces the field in the half space c < —2d. The field can/ be obtained from the 
free-space integral 1 ‘ 
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E(f)= ~Jh 


rin 

Jo Jv 


(r~r' - u(1 - 2d)) • J,•(f/,/)<W</^ , 


(63) 


Again, for a dipole source, the image current function can be written without t jie ( vari- 
able, because of the special 6(z' *'• A)- dependence of the original current source. Hence, 
to calculate the field, instead of both 1 and z' integrations we need only make the ^ in- 
tegration taking 4>[z' + 2d - A) instead of and integrating ~ from 1 A - 2d to -t~oc. 
From the transmission image functions we can see that the final image consists of dipoles 
and partially overlapping line sources alqng 'the positive z axis, each starting at the points 
A + (2m -i- 1 )2d, Thus, the transmission image current starts at. the height A-f 2d — 2d = A, 
or from the original source point. 

The images corresponding to magnetic sources can be obtained easily from the duality 
transform / — * 7 m ,c —» p. Note that c and p are relative quantities so that minus signs do 
not occur in this simple duality transform. 


ASYMPTOTIC TESTS - s4^/ 

To gain confidence in the above expressions, some of which have rather complicated 
appearance, it is necessary to perform asymptotic tests to see that correct limit cases are 
obtained. In fact, we should check that the image solution corresponds to 

(A) for c = 1 and p = 1, the free-space problem 

(B) for e — r oo, the perfect conductor plane problem 

(C) for p — y oo, the perfect magnetic conductor plane problem 

(D) for d — r 0, the free-space problem ~ ~ - 

(E) for d — * oo, the Sommerfeld half-space problem. 

The case (A) is easily seen to be satisfied, because B — >/fie - 1 = 0, whence the arguments 
of all Bessel functions in the definitions (21), ( 22) are zero, making the F m functions vanish. 
Thus, from (36), (39), (45), (50), (56) and (58) we see that all, the other images vanish 
except the delta function 6(t-2d) in the transmission image expressions, which corresponds 
to the original current as a transmission source and there is no reflection sources. 

The cases (B) and (C) are more complicated. In case (B), E = 1 and from (27) we can 
derive (he limit function needed here, 


c; A n (i) = (-i)% 0 - 


(64) 


Considering the VED source and inserting (64) in (34), it is seen that, the sum expressions 
all vanish, because Fo(0, /) = 0 and the terms in brackets cancel. Thus, only the first delta 
term survives giving f™ (t) = -6(1), which together with the mirror operation, c in (59), 
leaves us with the positive image of the vertical dipole. This corresponds to the mirror 
image due to a conducting plane. Also, it is readily seen from (38) that the transmission 
image of a VED vanishes. For the HED, the outcome is not so simply seen, because the 
image functions do not explicitly depend on e. The functions F m (c,t ), however, depend 
on the parameter B = y/fie - 1, which becomes large. Considering the image function 
expressions, we see that they can be written in the form Bf(Bi), and denoting a; = Bt , 
the field integral is of the form 

r G(i)Bf(Bt)di= r G(^)f(x)dx ^ G(0) r f(x)dx, (65) 

Jo Jo H Jo 
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as B -*•* pc. This means that the image can be replaced by a point source at t — 0, or 
: — ~ h , /What the point source is, is seen after the integral is developed. It is not difficult 
to evaluate these integrals if the identity (20) is invoked. In fact, writing for /3 = 0, /3j = B, 
we have j 


.1 

/Op 

/ (F w (e, 

h 


0 + 6 m o6(1 ~c))d1 ~ 


and taking a vanishingly small negative phase angle for B , the right-hand side vanishes, 
whence the integral of F m (c,t) for m > 0 is seen to vanish and for m = 0 give the value 
~*1. With this, checking all the components, we are able to show that for £ — ♦ oo, the 
image current source can be replaced by 


i f<X> 

Mr, t)dt = 

J 0 


This is' th<t mirror image due to the p erfectly conducting plane at z ==j 0. In the same 
manner, tht’ transmission image can be seen to vanish in this limit, Case j(C) follows from 
(B) through the duality transformation. - 

in case (D). tile functions Ff.(n2d,1) can be replaced by -whence in all image 

functions in the summations over,tlie Judex n most of the terms are cancelled. This makes 
the reflection image functions M (i), fo(0 vanish together with the transmission 
image function ho(t), leaving only — hT M (t) — 6(1 — 2 4 ), which corresponds to the 

primary source of the free-space problem. j ' ■■ 

Finally, in case (E), we should end up in the image of the Sommerfeld problem given 
in [1], To conform with the notation in [1], a change of the variable from t to p = jBi 
must be made. Then, writing! from (24), we have 


^.( 0,0 = = jb(-i) 




which gives us the following relation of the present notation and the function f a (p) of [1]: 


Y t Ci(E)F k (0,t) = -jBMp). (69) 

k=0 

In fact, the Sommerfeld image functions can be seen to result for the reflection image in 
the first interval from 1 = 0 to / = ~~4d for any value of d. For d — * oo, obviously other 
images are pushed to infinity and have no effect and the whole problem is reduced to the 
Sommerfeld problem. Concentrating on the interval t w 0... — 4 d, functions of the type 
6(1 — md) and /)’ with m > 0 can be discarded, since they are zero up to 1 = md, 

and only terms with 6(1) and Fi.(0,i) survive. Thus, we can write in this interval from 
(34), (69): 

; f™(t) = -jB (^y*(p) + Mp)) , (70) 

f TE (t)=jB^~~6( P )^U(p)y (71) 

"" - ~’ B + i < 72 > 

Because f' Q (1) - (jB)~g'(p), we have from the general image of (60) 
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i ^'(»V) - jBJj{r,p), or Ji(r,t)d1 = Jj(r t p)d.p, (73) 

where the left-hand side,* refer 1o the present theory and the right-hand sides to the cor- 
responding function pf the Sommerfeld problem as given in [1, p.1030]. Thus, the image 
function in the slab problem in the interval I = 0..i h 4<f is only dependent on the first 
interface, corresponding to the first reflected partial ray in Fig. 1. 

As an encore, we briefly show how the delta function terms in image functions can be 
obtained directly from the reflection coefficient expressions as a limit when /J — > oo. This 
corresponds to the simpler analysis of the delta function part of the VED image-given in 
[1, Part If], Studying (2) for this limit, whence /?] «-» f3, we see that because r™ — > — E, 
from the limit expression of R™ , eqn. (2), denoted by Rj M , we have the equation 

(1 - E 2 c~ 4j0d )Rj M = -J5( 1 - e~ 4j0d ) _._J : (74) 

This can be readily Laplace transformed to a difference equation for the corresponding 
image function i : 5 tl 


E 2 fJ M (t - Ad)V(i - Ad) =r\— 

The solution of this is obtained quite straightforwardly: 


() + E6(1 - Ad). 


oc 

■E6(1) + J7((-JS) 2n+1 - (~E ) 2n ~ *) 6(t - And), 


which is exactly the delta-function part of (34), 


GUIDED MODES 

It was seen in the analysis of the microstrip geometry [3], [4] that some of the image 
functions would not, converge along the real z axis and ’the reason was found to be in 
the nonradiating guided modes. After extracting the exponentially diverging image terms 
corresponding to these modes, the remaining image was convergent. The number of nonra- 
diating modes was dependent on the frequency and medium parameters, but there was at. 
least, one such mode ( TM mode) for the microstrip. The same applies for the slab except 
that there exist at least, two such modes (TE and' TM) in the slab. After extraction of 
these modes, the image functions are convergent. 

The modes can be classified in four groups in terms of TE/TM polarization and prop- 
erties of symmetry /antisymmetry of the transverse modal electric field. Antisymmetric 
modes correspond to those existing in the microstrip structure, whereas the symmetric 
ones exist in the microstrip with magnetic conducting ground plane. 

The modes are obtained from the poles of the reflection coefficients, all of which appear 
in the expression (12). The poles (3 = / 3 P , j3\ = f3 ip are obtained as roots of the four 
equations 

r TE <T M (i3,i3 l )e-j 0 ' 2d = ±1, ' (77) 

These equations can also be written as 


r TM c ~j0\ld _ or 

Jl-jlanUhd) = - A 

(77a) 


u>ee 0 

we 0 

r TM c ~j0\ 2d or 

A-\-jcot(0,d)) = -A 

'(776) 


cctfo 

wto 


~ h or 

7. 

(77c) 

;j ■■ r TE e~^ 2d =-l, or 

I . ; ' • ■ 


(77 d) 


The ri^lit-hand equations of (77a), (77c) correspond to transverse impedance interpretation 
of the problem with a conducting plane iiitlie middle of the slab and those of (77b), (77d) 
the same with a magnetically conducting plane. Thus, the problem is reduced to finding 
the modes in microstrip and magnetic microstrip structures, Each root of (77 a) ■ (77 d) 
corresponds to a guided mode in the slab. For sufficiently low values of Bd — y/fit — 1 led, 
only one TE and one TM mode have field dependencies that are decaying in the £ direction 
(imaginary 0) and can be considered as nonradiating waveguide inodes in the slab. These 
modes correspond to lowest solutions of (77a) and (77d), which can be obtained with the 
method described in [4j. Defining q ~ 0\ p d, we have for both the TM and TE wave J 

.. _ \ .. . \ 

■ v j : qtanq = yyJ(Bd) 2 - q 2 , (78) 

where 7 = t for the TM wave and 7 = p for the TE wave. The lowest solution for small 
values of Bd can easily be obtained, starting from go = Bd > n the iteration formula [4] 


0»+l = BdyJ 1 - ( q n 1anq n /~fBd ) 2 , (79) 

'! t 

after which 0p is obtained from (77a)-(77d). For real Bd, the graphical solution of Fig.2 
can also be applied. In Fig. 3, values of a p d — J0 p d for different Bd and e or p values are 
given, Corresponding to the respective cases of TM and TE waveguide modes. 

Because the guided mode is nonradiating, the corresponding pole satisfies lm[0 p \ < 0, 
which corresponds to an exponentially diverging reflection image current. This must be 
extracted from the image function to obtain a converging source. For this purpose, residue 
of the reflection coefficient at the pole must be calculated. We can write the TE and TM 
reflection coefficients in the same form 

with 0 P = ~ja p , the basic pole and j4p,|the corresponding residue j .. 



± 


; 7p?p. . 

I'B 2 - j0 P d(y 2 0 2 - 0\ p ) 


0ip ± y/ppA- B 2 . 



For the TE wave,! the uppersign and y — p must be taken, while for the TM wave, lower 
sign and 7 — e applies. 

In the limit Bd — > 0, the lowest root is g = 0i p d % Bd, whence 0 p d —ja p d = 
(—j0l p d/'y)tan0] p d —j(Bd) 2 / 7. In this case, the residue (81) reduces to A p % ±0 P as 
^fjB 2 d/y. The nondiverging part of the image function is /(/) — jApe^p 1 sg /(/) 7: 
(B 2 d/y)ixp(B 2 di/y). 
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The function jApe^, subtracted from the image function, diverges for real 1 and 
must he moved into another position on the complex t plane, as described in J3J and [4]. 


CALCULATION OF IMA'^E FUNCTIONS 


Finally, let us study numerical properties of the different image functions. The function 
F n (c, t) was defined in (21 ), (22) in terms of modified Bessel functions. For small argument 
values, the following power series, based on Taylor expansion of Bessel functions, can be 
applied [4j: _ • ' 


F 0 (c, i ) 


B 2 r 

2 


k k 

x y 


S *|(* + 1 )! 

B^* k { v 2n+*-l \ ,, 

F«lc,l) ~ "'2 S i! (( 2 n + * + !)! ~ (2n + * - l)t) ' '' ~ C *’ " f 0 


C('-c), 


with 


x = 


B{i + c) 


y- 


(2 n. + k 
B{t — c) 


(82) 


2 2 

Another method suitable for microcomputer use is to apply polynomial approximations 
for the modified Bessel functions /q(*), I\{x) and to compute functions with other indices 
from the recursive formula [7], after which (21), (22) can be applied. 

Examples of the normalized function jfF n (c,t) are depicted in Fig.4 for different val- 
ues o f n a nd Be, Because the modified Bessel functions J n ( a:) diverge asymptotically as 
e z / \/2irx for all n, the F n functions diverge for /->ixas 




(83) 


It was seen in [4] for the microstrip geometry that, in spite of the divergence of the indi- 
vidual F n terms in the image function series expression, the whole function may be well 
convergent wit hin the limit of computer accuracy. 

The functions Fm(c,t) can be grouped for better convergence. For example, denoting 
r m = By/t 2 — (4 rod)?, the special TM function of (36) can be written as 


— fTM 

B J 


(<) = - - B(t 

+ {l a (Bi)U(t)-(B(t~4d)) 


4d) h£il V{ri A 

n > 


zh( r \) 


) 


( B{i~m 


hh(7i) 


u(ri) - m - s d)) 


( 


h h(n) 
T$ 


+ (B(t - 


U(ri)-(B(1 


8 d)) 7 ^W 2 


) 


(84) 
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For values of Bd satisfying Bd < ir/2 there only exists one nonradiating guided TM 
inode, whose parameters can be obtained from (78). For small values of Bd , we can 
write the two-term approximation of the function (84) with the pole term removed in the 
following simple form; 


B 


(,f™(t) - jApc’W) 


B tI t , 1 IB 


(f) : 


(( -4d) 2 {t + 8d) 


\U(t)~U(i~4d)\ + ^(2 BdfU(t). (85) 


This shows us clearly how the image function tends to a simple triangular pulse within the 
interval (0,4d) and is small elsewhere, 

Values from (84) both with and without the guided mode term were calculated with a 
PC using simple routines for the modified Bessel functions [9], Results are given in Figs. 
5 and 6. The convergence is seen to correspond to the expression (85) for the dashed 
function despite the rapid divergence (83) of tne F„ functions. Thus, electrically thin slabs 
can be analyzed extremely easily using the present theory. 

Slightly more efTort is needed for the evaluation of the general image functions like 
(34). We calculated values for the image function for e = 2 for two normalized 

thicknesses of the slab, 4 Bd = 1 and 4 Bd. — 2, Figs. 7,8. Here again, in spite of the inherent 
divergence of the modified Bessel functions, the image function without the guided mode 
term (dashed lines) was seen to converge for values of t small enough where the Bessel 
function routines were accurate. 

In addition to the delta function terms, there are also step discontinuities at the values 
t = 4md of the argument. It is not difficult to find exact expressions for these step 
discontinuities of the image function, in fact, from the k ~ 0 terms of (34) we obtain 


Uf™w - 


E 1M '(I - E?)Bimd at i - 4md. 


(86) 


There are obviously no step discontinuities when e == 1 (Figs. 5, 6) ort = oo. The largest 
step will occur at m = (c + l) 2 /4 c. The first step is largest for t < 5.8. With the present 
PC and simple Bessel function routines, the applicability of the expressions, seems to end 
at about e % 4 with moderate thicknesses of the slab (only first imodes propagating). Thus, 
the method seems to be best suited at the moment for microstrip antenna structures, where 
low ( is desired. , - — r j 

The function f TE (1) is of the same form as the corresponding f™ (1) function and 
can be obtained from it aftep the interchange t p and; change of sign. Evaluation of 
the more complicated /o(f ) function and the corresponding transmission image functions 
must be left to a future work. It is hoped that an easier analytic form for the image 
functions could be found with the modified Bessel functions replaced by functions with 
intrinsic convergence properties.; Since the exact expressions ! have been already derived, 
it seems to be a question of finding suitable analytic transformations for these functions. 
The residue series method applied in [3] is, of course, one such possibility, leading to a 
series of exponent ial functions, but it has the inconvenience of requiring a determination 
of a set of poles from solutions of a transcendental equation and, thus, does not yield the 
luxury of an explicit expression. 
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CONCLUSION 


The exact image met hod has been extended to problems involving a slab with ( and/or 
fi different frtmi those of the surrounding space, Both reflection and transmission image 
expressions hive been derived and their validity studied with asymptotic tests. Numerical 
convergence has been demonstrated with calculations, With jibe present experience, using 
a PC and simple routines for the modified Bessel functions, the method appears to work 
for values of c < 3 and moderate thicknesses of the slab and, thus, be of interest for the 
analysis of microstrip antennas. Further effort should be directed to transforming the 
present exact expressions into other forms for the extension of their numerical range of 
validity in microcomputer calculations. 
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FIGURE C APTIONS 

■ ! 

Figure I, Schematic representation of the partial waves in the slab problem showing 
the multiple reflections and transmissions resulting in the total reflection and transmission 
coefficients for a plane wavt* (Fourier space component) arising from the original source. 
The rays, act ually parallel to t lie - axis, are drawn at an angle for clarity. 

Figure 2. Graphical method foi‘ the construction of roots corresponding to the lowest 
TE and TM guided modes in the slab of thickness 2d. 

Figure 3. Values of the attenuation factor of a guided wave in the slab, a p d for different 
values of the parameters Bd. apd 7. Both TE (with 7 = p) and TM (with 7 — c) cases 
are covered by the same diagram. 


Figure 4. Normalized function F n (c,i) for different values of n and Be, 

Figure 5. Normalized image function f™(t) (solid line), as calculated from (84), and 
the same with the residue (guided mode) image term extracted (dashed line) in the case 
f = 1, fi 7^ 1, 4 Bd ~ 1. The corresponding function for fi = 1, e ^ 1, 4 Bd = 1 is 

obtained by merely changing the sign. 

Figure 6. The same as Fig. 5, with ABd = 2. 

-Figure 7. Normalized image function without the delta function terms, calcu- 

lated from (34a) (solid line) and the same with the residue (guided mode) term extracted 
(dashed line), for e *= 2, ABd = 1. The corresponding function for fi = 2, 4 Bd = 1 

is obtained by merely changing the sign. 

Figure 8. The same as Fig. 7, with ABd — 2. 
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